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Search for massive supersymmetric particles decaying to many jets using the ATLAS 

detector in pp collisions at y / s = 8 TeV 

The ATLAS Collaboration 

Results of a search for decays of massive particles to fully hadronic final states are presented. 

This search uses 20.3 fb -1 of data collected by the ATLAS detector in yfs = 8 TeV proton-proton 
collisions at the LHC. Signatures based on high jet multiplicities without requirements on the 
missing transverse momentum are used to search for Jt-parity-violating supersymmetric gluino pair 
production with subsequent decays to quarks. The analysis is performed using a requirement on 
the number of jets, in combination with separate requirements on the number of b-tagged jets, as 
well as a topological observable formed from the scalar sum of the mass values of large-radius jets 
in the event. Results are interpreted in the context of all possible branching ratios of direct gluino 
decays to various quark flavors. No significant deviation is observed from the expected Standard 
Model backgrounds estimated using jet-counting as well as data-driven templates of the total-jet- 
mass spectra. Gluino pair decays to ten or more quarks via intermediate neutralinos are excluded 
for a gluino with mass m,g < 1 TeV for a neutralino mass m^o = 500 GeV. Direct gluino decays to 
six quarks are excluded for < 917 GeV for light-flavor final states, and results for various flavor 


hypotheses are presented. 

PACS numbers: 12.60.Jv,12.38.Qk,11.30.Pb,13.87.-a 

I. INTRODUCTION 

Supersymmetry (SUSY) [1-9] is a theoretical exten¬ 
sion of the Standard Model (SM) which fundamentally 
relates fermions and bosons. It is an alluring theoret¬ 
ical possibility given its potential to solve the natural¬ 
ness problem [10-15] and to provide a dark-matter can¬ 
didate [16, 17]. Partially as a result of the latter pos¬ 
sibility, most searches for SUSY focus on scenarios such 
as a minimal supersymmetric standard model (MSSM) 
in which .R-parity is conserved (RPC) [18-21], In these 
models, SUSY particles must be produced in pairs and 
must decay to a stable lightest supersymmetric particle 
(LSP). With strong constraints now placed on standard 
RPC SUSY scenarios by the experiments at the Large 
Hadron Collider (LHC), it is important to expand the 
scope of the SUSY search program and explore models 
where R-parity may be violated and the LSP may decay 
to SM particles, particularly as these variations can alle¬ 
viate to some degree the fine-tuning many SUSY models 
currently exhibit [22]. 

In i?-parity-violating (RPV) scenarios, many of the 
constraints placed on the MSSM in terms of the allowed 
parameter space of gluino ( g ) and squark ( q ) masses 
are relaxed. The reduced sensitivity of standard SUSY 
searches to RPV scenarios is due primarily to the high 
missing transverse momentum (A™ lss ) requirements used 
in the event selection common to many of those searches. 
This choice is motivated by the assumed presence of two 
weakly interacting and therefore undetected LSPs. Con¬ 
sequently, the primary challenge in searches for RPV 
SUSY final states is to identify suitable substitutes for 
the canonical large A™ lss signature of RPC SUSY used 
to distinguish signals from background processes. Com¬ 
mon signatures used for RPV searches include resonant 
lepton pair production [23-25], exotic decays of long-lived 


particles, and displaced vertices [26-29]. 
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Figure 1. Diagrams for the benchmark processes considered 
for this analysis. The solid black lines represent Standard 
Model particles, the solid red lines represent SUSY partners, 
the gray shaded circles represent effective vertices that in¬ 
clude off-shell propagators (e.g. heavy squarks coupling to a 
Xi neutralino and a quark), and the blue shaded circles rep¬ 
resent effective RPV vertices allowed by the baryon-number- 
violating A" couplings with off-shell propagators (e.g. heavy 
squarks coupling to two quarks). 

New analyses that do not rely on E™ lss are required 
in order to search for fully hadronic final states involving 
RPV gluino decays directly to quarks or via Xi neutrali¬ 
nos as shown in the diagrams in Fig. 1. Cases in which 
pair-produced massive new particles decay directly to a 
total of six quarks, as well as cascade decays with at least 
ten quarks, are considered. Three-body decays of the 
type shown in Fig. 1 are given by effective RPV vertices 
allowed by the baryon-number-violating A " couplings as 
described in Sec. II with off-shell squark propagators. 
This analysis is an extension of the search conducted at 
y/s = 7 TeV for the pair production of massive gluinos, 
each decaying directly into three quarks [30]. 

The diagrams shown in Fig. 1 represent the bench¬ 
mark processes used in the optimization and design of 
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the search presented in this paper. The extension to con¬ 
sidering cascade decays of massive particles creates the 
potential for significantly higher hadronic final-state mul¬ 
tiplicities and motivates a shift in technique with respect 
to previous searches. Therefore, the analysis is extended 
to look for events characterized by much higher recon¬ 
structed jet multiplicities as well as with event topolo¬ 
gies representative of these complex final states. Two 
complementary search strategies are thus adopted: a jet¬ 
counting analysis that searches for an excess of >6-jet or 
>7-jet events, and a data-driven template-based analy¬ 
sis that uses a topological observable called the total-jet- 
mass of large-radius (large- R) jets. The former exploits 
the predictable scaling of the number of n-jet events 
(n = 6, 7) as a function of the transverse momentum 
(pt) requirement placed on the n th leading jet in px for 
background processes. This analysis is sensitive to the 
models presented here because this scaling relation differs 
significantly between the signal and the background. The 
latter analysis uses templates of the event-level observ¬ 
able formed by the scalar sum of the four leading large -R 
jet masses in the event, which is significantly larger for 
the signal than for the SM backgrounds. 

This paper is organized as follows: Sec. II describes the 
motivation and theoretical underpinnings of the bench¬ 
mark processes used in this analysis. Section III and 
Sec. IV present details of the detector, the data collection 
and selection procedures, and the Monte Carlo (MC) sim¬ 
ulation samples used for this search. The physics object 
definitions used to identify and discriminate between sig¬ 
nal and background are described in Sec. V. The details 
of the methods are separated for the two analyses em¬ 
ployed. The jet-counting analysis is presented in Sec. VI, 
while the total-jet-mass analysis using more advanced ob¬ 
servables is presented in Sec. VII. The combined results 
of this search and the final sensitivity to the benchmark 
processes are then described in Sec. VIII. The results 
using the total-jet-mass analysis are presented first, in 
Sec. VIII A, as they only apply to the ten-quark final 
states. The jet-counting analysis additionally yields in¬ 
terpretations across the flavor structure allowed by the A" 
couplings. This comprehensive set of results is presented 
in Sec. VIIIB. Comparisons between the two analyses are 
then made in Sec. VIIIC. 


II. R-PARITY-VIOLATING SUPERSYMMETRY 
AND BARYON-NUMBER VIOLATION 

The benchmark model used to interpret the results of 
the search for high multiplicity hadronic final states is the 
baryon-number-violating RPV SUSY scenario. The RPV 
component of the generic supersymmetry superpotential 
can be written as [31, 32]: 

W K — 2 ^ijkLiLj Ek T A- L, Cjj /-A + 

I KjkUiDjDk+mUHz, (1) 

where i,j,k = 1,2,3 are generation indices. The gen¬ 


eration indices are sometimes omitted in the discussions 
that follow if the statement being made is not specific 
to any generation. The first three terms in Eq. (1) are 
often referred to as the trilinear couplings, whereas the 
last term is bilinear. The Li, Qi represent the lepton 
and quark SU{ 2 )l doublet superfields, whereas Hi is the 
Higgs superfield. The Ej, Dj, and Uj are the charged lep¬ 
ton, down-type quark, and up-type quark SU( 2)l singlet 
superfields, respectively. The Yukawa couplings for each 
term are given by A, A', and A", and n is a dimensionful 
mass parameter. In general, the particle content of the 
RPV MSSM is identical to that of the RPC MSSM but 
with the additional interactions given by Wft,,. 

Generically, the addition of Wjt p into the overall SUSY 
superpotential allows for the possibility of rapid pro¬ 
ton decay. The simultaneous presence of lepton-number- 
violating (e.g. A' y^ 0) and baryon-number-violating op¬ 
erators (A" y£ 0) leads to proton decay rates larger than 
allowed by the experimental limit on the proton lifetime 
unless, for example, [33] 
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where rriq is the typical squark mass. As a result, even 
when considering this more generic form of the SUSY 
superpotential by including Wft , it is still necessary to 
impose an ad hoc, albeit experimentally motivated, sym¬ 
metry to protect the proton from decay. It is generally 
necessary that at least one of A, A', A" be exactly equal to 
zero. Consequently, it is common to consider each term 
in Eq. (1) independently. In the case of nonzero A and A 7 , 
the typical signature involves leptons in the final state. 
However, for A" k y^ 0, the final state is characterized by 
jets, either from direct gluino decay or from the cascade 
decay of the gluino to the lightest neutralino (Xi), as also 
considered here. Because of the structure of Eq. (1), sce¬ 
narios in which only \"- k y^ 0 are often referred to as 
UDD scenarios. 

Current indirect experimental constraints [34] on the 
sizes of each of the UDD couplings A"- fc from sources other 
than proton decay are valid primarily for low squark 
masses, as suggested by Eq. (2). Those limits are driven 
by double nucleon decay [35] (for A" 12 ), neutron oscilla¬ 
tions [36] (for A] 7 13 ), and Z boson branching ratios [37]. 

Hadron collider searches are hindered in the search for 
an all-hadronic decay of new particles by the fact that the 
SM background from multi-jet production is very high. 
Nonetheless, searches have been carried out by several 
collider experiments. The CDF Collaboration [38] ex¬ 
cluded gluino masses up to 240 GeV for light-flavor mod¬ 
els. The CMS Collaboration [39] excludes such gluinos 
up to a mass of 650 GeV and additionally sets limits on 
some heavy-flavor UDD models. The ATLAS Collabo¬ 
ration [30] has also previously set limits in a search for 
anomalous six-quark production, excluding gluino masses 
up to 666 GeV for light-flavor models. The search pre¬ 
sented here uniquely probes the flavor structure of the 
UDD couplings and employs new techniques both in anal- 
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ysis and theoretical interpretation. 

III. THE ATLAS DETECTOR 

The ATLAS detector [40] provides nearly full solid 
angle coverage around the collision point with an in¬ 
ner tracking system covering the pseudorapidity 1 range 
\r/\ < 2.5, electromagnetic (EM) and hadronic calorime¬ 
ters covering |?;| < 4.9, and a muon spectrometer covering 
|r;| <2.7. 

The ATLAS tracking system is composed of a silicon 
pixel tracker closest to the beamline, a microstrip silicon 
tracker, and a straw-tube transition radiation tracker. 
These systems are layered radially around each other 
in the central region. A thin solenoid surrounding the 
tracker provides an axial 2 T field enabling measurement 
of charged-particle momenta. 

The calorimeter, which spans the pseudorapidity range 
up to h| = 4.9, is comprised of multiple subdetectors 
with different designs. The high granularity liquid argon 
(LAr) electromagnetic calorimeter system includes sepa¬ 
rate barrel (hi < 1-475), endcap (1.375 < hi < 3.2), and 
forward subsystems (3.1 < hi < 4.9). The tile hadronic 
calorimeter (|? 7 | < 1.7) is composed of scintillator tiles 
and iron absorbers. As described below, jets used in the 
analyses presented here are typically required to have 
hi < 2.8 such that they are fully contained within the 
barrel and endcap calorimeter systems. 

A three-level trigger system is used to select events 
to record for offline analysis. The level-1 trigger is im¬ 
plemented in hardware and uses a subset of detector in¬ 
formation to reduce the event rate to a design value of 
at most 75 kHz during 2012. This is followed by two 
software-based triggers, level-2 and the event filter (col¬ 
lectively called the high-level trigger), which together re¬ 
duce the event rate to a few hundred Hz. The primary 
triggers used in this analysis collected the full integrated 
luminosity of the 8 TeV dataset with good efficiency for 
the event selections described in this paper. 

IV. DATA AND MONTE CARLO SAMPLES 

The data used in this analysis correspond to 20.3 ± 
0.6 fb _1 [41, 42] of integrated luminosity taken during pe¬ 
riods in which the data satisfied baseline quality criteria. 
Further details of the event selections applied, including 
the ATLAS data quality criteria and trigger strategy, are 


1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the center of the detector 
and the 2 -axis along the beam pipe. The x-axis points from the IP 
to the center of the LHC ring, and the y -axis points upward. Cylin¬ 
drical coordinates (r. <!>) are used in the transverse plane, 0 being 
the azimuthal angle around the beam pipe. The pseudorapidity is 
defined in terms of the polar angle 8 as y = — lntan(0/2). 


given in Sec. V. The primary systems of interest in these 
studies are the electromagnetic and hadronic calorimeters 
and the inner tracking detector. The data were collected 
with triggers based on either single-jet or multi-jet sig¬ 
natures. The single-jet trigger selection has a transverse 
momentum threshold of 360 GeV using a large -R anti-fc t 
jet definition [43] with a nominal radius of R = 1.0 within 
the high-level jet trigger. The multi-jet trigger selection 
requires at least six anti-fc* R = 0.4 jets with a nominal 
Pt threshold of 45 GeV in the high-level trigger. Data 
collected using several additional multi-jet requirements 
(from three to five jets) are also used for background es¬ 
timation studies. 

Multiple simultaneous proton-proton (pp) interac¬ 
tions, or pile-up, occur in each bunch crossing at the 
LHC. The additional collisions occurring in the same and 
neighboring bunch crossings with respect to the event of 
interest are referred to as in-time and out-of-time pile- 
up, respectively, and are uncorrelated with the hard- 
scattering process. 

The benchmark RPV SUSY signal processes of both 
the six-quark and ten-quark models (see Sec. I) were sim¬ 
ulated using HERWIG++ 6.520 [44] for several gluino and 
neutralino mass hypotheses using the parton distribution 
function (PDF) set CTEQ6L1 [45, 46]. For both models, 
all squark masses are set to 5 TeV and thus gluinos decay 
directly to three quarks or to two quarks and a neutralino 
through standard RPC couplings. In the ten-quark cas¬ 
cade decay model, the neutralinos each decay to three 
quarks via an off-shell squark and the RPV UDD decay 
vertex with coupling A". In this model, the neutralino is 
the lightest supersymmetric particle. 

Samples are produced covering a wide range of both 
Trig and to - o. In the six-quark direct gluino decay model, 
the gluino mass is varied from 500 to 1200 GeV. In the 
case of the cascade decays, for each gluino mass (400 GeV 
to 1.4 TeV), separate samples are generated with multi¬ 
ple neutralino masses ranging from 50 GeV to 1.3 TeV. 
In each case, to-o < rrig. In order to ensure the re¬ 
sult has minimal sensitivity to the effects of initial state 
radiation (ISR), which could be poorly modeled in the 
signal samples, 2 the region with ( rrig — to- o) < 100 GeV 
is not considered. Due to the potentially large theoreti¬ 
cal uncertainty on the non-SM colorflow given by UDD 
couplings, results are presented for a single model of ra¬ 
diation and no systematic uncertainty is assigned for this 
effect, further justifying the unevaluated region described 
above. All possible A"- fe flavor combinations given by the 
structure of Eq. (1) are allowed to proceed with equal 
probability. As discussed in Sec. VIII, the analysis main¬ 
tains approximately equal sensitivity to all flavor modes. 
All samples are produced assuming that the gluino and 


2 HERWIG-| — , which is used for signal simulation, is not ex¬ 
pected to model additional energetic jets from ISR well because the 
leading-order evaluation of the matrix element is only performed 
for the 2 —» 2 particle scattering process. 
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neutralino widths are narrow and that their decays are 
prompt. Cross-section calculations are performed at 
next-to-leading order in the strong coupling constant, 
adding the resummation of soft gluon emission at next- 
to-leading-logarithmic accuracy (NLO+NLL) [47-51]. 

Dijet and multi-jet events, as well as top quark pair 
production processes, were simulated in order to study 
the SM contributions and background estimation tech¬ 
niques. In the case of the vastly dominant background 
from SM jet production, several MC simulations were 
compared with data for the suitability of their descrip¬ 
tions of jet and multi-jet kinematic observables and 
topologies. For signal region selections that use ^-tagging 
(the identification of jets containing R-hadrons), other 
backgrounds such as tt, single top, and W/Z +jets be¬ 
come significant as well. These other backgrounds are 
estimated directly from the simulation. 

In order to develop the data-driven background esti¬ 
mation techniques for multi-jet events from QCD pro¬ 
cesses, comparisons are made among various generators 
and tunes. In the case of the jet-counting analysis, the 
ATLAS tune AUET2B LO** [52] of PYTHIA 6.426 [53] 
is used in estimating the rate of n-jet events (where 
n = 6, 7) as a function of the jet px requirement on the 
n th jet. For the total-jet-mass analysis, SHERPA 1.4.0 [54] 
is used to develop and test the method. For the SHERPA 
multi-jet samples, up to three partons are included in 
the matrix-element calculation and no electroweak pro¬ 
cesses are included. Heavy (c and b ) quarks are treated 
as massive. The next largest background after multi-jets 
is fully hadronic tt production, which is also simulated 
with SHERPA 1.4.0 and is used to estimate any background 
contamination in the control and signal regions defined 
in the analysis. 

The jet-counting and total-jet-mass analyses use dif¬ 
ferent multi-jet generators because of the different ap¬ 
proaches to the background estimation employed by each 
analysis. The low-to-high jet-multiplicity extrapolation 
of the jet-counting analysis, described in Sec. VIA, fa¬ 
vors a generator that treats the production of an addi¬ 
tional jet in a consistent manner, such as PYTHIA, rather 
than a generator that treats the multileg matrix element 
separately from the additional radiation given by a sepa¬ 
rate parton shower model. In contrast, the total-jet-mass 
analysis uses the multi-jet simulation only to test the 
background estimation method and optimize the anal¬ 
ysis as described in Sec. VIIA and Sec. VIIB, and uses 
SHERPA as it provides a better description of jet substruc¬ 
ture variables, such as the jet mass used in this analysis. 

The ATLAS simulation framework [55] is used to pro¬ 
cess both the signal and background events, including a 
full GEANT4 [56] description of the detector system. The 
simulation includes the effect of both in-time and out-of- 
time pile-up and is weighted to reproduce the observed 
distribution of the average number of collisions per bunch 
crossing in the data. 


V. PHYSICS OBJECTS AND EVENT 
PRESELECTION 

A. Data quality criteria 

The data are required to have met criteria designed to 
reject events with significant contamination from detec¬ 
tor noise, noncollision beam backgrounds, cosmic rays, 
and other spurious effects. The selection related to these 
quality criteria is based upon individual assessments for 
each subdetector, usually separated into barrel, forward 
and endcap regions, as well as for the trigger and for each 
type of reconstructed physics object (i.e. jets). 

To reject noncollision beam backgrounds and cosmic 
rays, events are required to contain a primary vertex 
consistent with the LHC beamspot, reconstructed from 
at least two tracks with transverse momenta plf ack > 
400 MeV. Jet-specific requirements are also applied. All 
jets reconstructed with the anti -k t algorithm using a ra¬ 
dius parameter of R = 0.4 and a measured > 20 GeV 
are required to satisfy the “looser” requirements discussed 
in detail in Ref. [57]. This selection requires that jets de¬ 
posit at least 5% of their measured total energy in the 
electromagnetic (EM) calorimeter as well as no more than 
99% of their energy in a single calorimeter layer. 

The above quality criteria selections for jets are ex¬ 
tended to prevent contamination from detector noise 
through several detector-region-specific requirements. 
Jets with spurious energy deposits in the forward 
hadronic endcap calorimeter are rejected and jets in the 
central region (|r;| < 2.0) that are at least 95% contained 
within the EM calorimeter are required not to exhibit any 
electronic pulse shape anomalies [58]. Any event with a 
jet that fails the above requirements is removed from the 
analysis. 

B. Object definitions 

Jets are reconstructed using the anti-fct algorithm with 
radius parameters of both R = 0.4 and R = 1.0. The 
former are referred to as standard jets and the latter as 
large-R jets. The inputs to the jet reconstruction are 
three-dimensional topological clusters [59]. This method 
first clusters together topologically connected calorimeter 
cells and classifies these clusters as either electromagnetic 
or hadronic. The classification uses a local cluster weight¬ 
ing (LCW) calibration scheme based on cell-energy den¬ 
sity and longitudinal depth within the calorimeter [60]. 
Based on this classification, energy corrections derived 
from single-pion MC simulations are applied. Dedicated 
corrections are derived for the effects of noncompensa¬ 
tion, signal losses due to noise-suppression threshold ef¬ 
fects, and energy lost in noninstrumented regions. An 
additional jet energy calibration is derived from MC sim¬ 
ulation as a correction relating the calorimeter response 
to the true jet energy. In order to determine these correc¬ 
tions, the identical jet definition used in the reconstruc- 
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tion is applied to particles with lifetimes greater than 
10 ps output by MC generators, excluding muons and 
neutrinos. Finally, the standard jets are further cali¬ 
brated with additional correction factors derived in-situ 
from a combination of 7 +jet, Z+ jet, and dijet balance 
methods [60]. 

No explicit veto is applied to events with leptons or 
U™ lss . This renders the analysis as inclusive as possible 
and leaves open the possibility for additional interpre¬ 
tations of the results. There is no explicit requirement 
removing identified leptons from the jets considered in an 
event. Calorimeter deposits from leptons may be consid¬ 
ered as jets in this analysis given that the data quality 
criteria described in Sec. V A are satisfied. A further con¬ 
sequence of these requirements is that events containing 
hard isolated photons, which are not separately identi¬ 
fied and distinguished from jets, have a high probability 
of failing to satisfy the signal event selection criteria. For 
the signals considered, typically 1 % of events fail these 
quality requirements. 

The standard jet px requirement is always chosen to 
be at least 60 GeV in order to reside in the fully efficient 
region of the multi-jet trigger. For the jet-counting anal¬ 
ysis selection (Sec. VI), a requirement of p^' > 80 GeV 
is imposed for each jet in most of the background control 
regions, and a higher requirement is used for the major¬ 
ity of the signal regions of the analysis. All jets used in 
this analysis are required to have |r;| < 2.8. The effect of 
pile-up on jets is negligible for the kinematic range con¬ 
sidered, and no selection to reduce pile-up sensitivity is 
included. 

In order to constrain specific UDD couplings to heavy 
flavor quarks, 6 -tagging requirements are also applied to 
some signal regions. In these cases, one or two standard 
jets are required to satisfy 6 -tagging criteria based on 
track transverse impact parameters and secondary vertex 
identification [61]. In simulated tt events, this algorithm 
yields a 70% (20%) tagging efficiency for real 6 -(c-)jets 
and an efficiency of 0.7% for selecting light quark and 
gluon jets. The 6 -tagging efficiency and misidentification 
are corrected by scale factors derived in data [61]. These 
jets are additionally required to lie within the range |? 7 | < 
2.5. 

The topological selection based on the total mass 
of large-I? jets (Sec. VII) employs the trimming algo¬ 
rithm [62]. This algorithm takes advantage of the fact 
that contamination from the underlying event and pile- 
up in the reconstructed jet is often much softer than the 
outgoing partons from the hard-scatter. The ratio of the 
px of small subjets (jets composed of the constituents of 
the original jet) to that of the jet is used as a selection 
criterion. The procedure uses a fc t algorithm [63, 64] to 
create subjets with a radius I? su b = 0.3. Any subjets 
with pxi/Px* < /cut are removed, where pxi is the trans¬ 
verse momentum of the i th subjet, and / cut = 0.05 is 
determined to be an optimal setting [65] . The remaining 
constituents form the trimmed jet, and the mass of the 
jet is the invariant mass of the remaining subjets (which 


in turn is the invariant mass of the massless topological 
clusters that compose the subjet). Using these trimming 
parameters, the full mass spectrum is insensitive to pile- 
up. 

The total-jet-mass analysis uses a sample from the 
high-py* single-jet triggers. A requirement that the lead¬ 
ing large-I? jet have pl^ > 500 GeV is applied to ensure 
that these triggers are fully efficient. 


VI. JET-COUNTING ANALYSIS 


A. Method and techniques 


The jet-counting analysis searches for an excess of 
events with > 6 or > 7 high px jets jets (with at least 
80 GeV), with > 0, > 1, or > 2 6 -jet requirements added 
to enhance the sensitivity to couplings that favor decays 
to heavy-flavor quarks. The number of jets, the px re¬ 
quirement that is used to select jets, and the number 
of 6 -tagged jets are optimized separately for each signal 
model taking into account experimental and theoretical 
uncertainties. 

The background yield in each signal region is estimated 
by starting with a signal-depleted control region in data 
and extrapolating its yield into the signal region using 
a factor that is determined from a multi-jet simulation, 
with corrections applied to account for additional minor 
background processes. This can be expressed as: 


N , 


_ (jy data jyMC 


n- jet — l/ v m -jet 


m- jet, Other BGs/ 

. n mc 

iy n- jet, Other BGs 


( N “St A 


(3) 


where the number of predicted background events with 
n jets (N n . j e t) is determined starting from the number 
of events in the data with m jets (A^A^). The extrap- 

N M 9 

olation factor, A'c* 1 , is determined from multi-jet simu- 

iV m-jet 

lation and validated in the data. This procedure is per¬ 
formed in exclusive bins of jet multiplicity. Since the sim¬ 
ulation is not guaranteed to predict this scaling perfectly, 
cross-checks in the data and a data-driven determination 
of systematic uncertainties are performed as described in 
Sec. VIC. It is assumed here that the simulation used 
for this extrapolation given by PYTHIA 6.426 predicts the 
relative rate of events with one additional order in the 
strong coupling constant in a consistent way across jet- 
multiplicity regions. This assumption comes from the 
behavior of the parton shower model used by PYTHIA to 
obtain configurations with more than two partons and is 
shown to be consistent with data in the measurement of 
multi-jet cross-sections [ 66 ]. Other models were studied 
and are discussed in Sec. VIC. 

Small corrections from other backgrounds (tt, single 
top, and W/Z+iet events) are applied based on estimates 
from the simulation. Without 6 -tagging, the contribution 
of events from these other backgrounds is less than 1 %. 
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Including two 5-tagged jets increases this relative contri¬ 
bution to as much as 10%. 


B. Signal and control region definitions 

Control regions are defined with m < 5, for which 
the background contribution is much larger than the ex¬ 
pected signal contributions from the benchmark signal 
processes. Extrapolation factors with n,m < 5 are used 
to validate the background model and to assign system¬ 
atic uncertainties. For n > 5, the expected signal con¬ 
tributions can become significant and an optimization is 
performed to choose the best signal region definitions for 
a given model. Signal regions are chosen with simulta¬ 
neous optimizations of the jet-multiplicity requirement 
(> 6 or > 7 jets), the associated transverse momentum 
requirement (80-220 GeV in 20 GeV steps), and the min¬ 
imum number of 5-tagged jets (> 0, > 1, or > 2) for 
a total of 48 possible signal regions. Alternative control 
regions are constructed from some n > 5 regions when 
the signal significance is expected to be low as described 
in Sec. VIC. Such regions are then excluded from the list 
of allowed signal regions. For a given signal model, the 
signal region deemed most effective by this optimization 
procedure is used for the final interpretations. The sig¬ 
nal regions chosen by the optimization procedure tend to 
pick regions with signal acceptances as low as 0.5% and 
as high as roughly 20%. 

Although other choices are also studied to determine 
background yield systematic uncertainties from the data, 
the background contributions are estimated in the fi¬ 
nal signal regions using extrapolations across two jet- 
multiplicity bins (n = m + 2). This choice leads to negli¬ 
gible signal contamination in the control regions used for 
this nominal prediction. 


C. Validation and systematic uncertainties 

Since the 3-, 4-, and 5-jet-multiplicity bins have mini¬ 
mal expected signal contamination they are used to vali¬ 
date the background model based on the MC simulation. 
The initial validation of the background prediction is per¬ 
formed by extrapolating the background from either the 
m = 3 or m = 4 jets control region into the n = 5 jets 
control region and comparing with the data. This com¬ 
parison is presented in Fig. 2, which shows the number 
of events passing a given jet px requirement with a 5- 
jet requirement. This procedure is shown to be accurate 
in the extrapolations to the 5-jet bin in data, both with 
and without the requirement of 5-tagging. The conclu¬ 
sion of this validation study is that Eq. (3) can be used 
with no correction factors, but a systematic uncertainty 
on the method is assigned to account for the discrepan¬ 
cies between data and the prediction in the control re¬ 
gions. This systematic uncertainty is assigned to cover, 
per p;p l bin, the largest discrepancy that is observed be¬ 


tween data and the prediction when extrapolating from 
either the 3-jet or 4-jet bins into the 5-jet control region, 
as well as from extrapolations to higher jet multiplicity 
as discussed below. 

Alternative MC models of extra-jet production such 
as those given by SHERPA, HERWIG++, and additional pa¬ 
rameter tunes in PYTHIA were studied and either did not 
satisfy the criterion that the model be consistent through 
control and signal regions (e.g. the model must not de¬ 
scribe the control regions with a matrix element calcula¬ 
tion and the signal regions with a parton shower model 
giving unreliable projections) or disagreed significantly 
with the data in the validations presented here. The in¬ 
ternal spread of predictions given by each of these back¬ 
ground models in various extrapolations is considered 
when assigning systematic uncertainties. In all cases, this 
spread is consistent with the systematic uncertainties ob¬ 
tained using PYTHIA in the manner described above. 

In addition to the extrapolation factor described by 
Eq. (3), it is possible to also study the extrapolation 
along the jet px degree of freedom. In this case, the n-jet 
event yield for a given high jet-px selection is predicted 
using extrapolation factors from lower jet-px selections 
determined from MC simulation. This method is tested 
exclusively in a low n-jet region for the high jet-px re¬ 
quirement and the spread is compared to the baseline 
systematic uncertainty, which is increased in case of dis¬ 
agreement larger than this baseline. 

Additional control regions can be constructed from ex¬ 
clusive 6-jet regions with low jet-px requirements. Any 
region with an expected signal contribution less than 10% 
for the nig = 600 GeV six-quark model is used as addi¬ 
tional control region in the evaluation of the background 
systematic uncertainties. These regions are used to en¬ 
sure that the jet-multiplicity extrapolation continues to 
accurately predict the event rate at higher jet multiplic¬ 
ities, as shown in Fig. 3(a), without looking directly at 
possible signal regions. This procedure allows the exclu¬ 
sive 6-jet, low jet-px region to be probed and shows that 
the jet-multiplicity extrapolations continue to provide ac¬ 
curate predictions at higher jet multiplicities. 

To extend this validation, a requirement that the av¬ 
erage jet pseudorapidity (|p|) > 1.0 is applied to create 
a high-pseudorapidity control region to reduce the signal 
contribution to a level of less than approximately 10% 
while retaining a reasonable number of events. Results 
of these extrapolations are shown for the exclusive 7-jet 
bin in Fig. 3(b). The largest deviations from the expected 
values are found to be a few percent larger than for the 
5-jet extrapolations. 

The uncertainty due to any mismodeling of contri¬ 
butions from backgrounds such as ft, single top, and 
W/Z+]et processes is expected to be small and is cov¬ 
ered by the procedure above since these contributions 
are included in the extrapolation. Therefore, any mis¬ 
modeling of these sources results in increased systematic 
uncertainty on the entire background model in this pro¬ 
cedure. 
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(a) > 0 b-tagged jets required (b) > 1 b-tagged jets required 



60 80 100 120 140 160 180 200 220 240 

Jet p T Requirement [GeV] 


(c) > 2 b-tagged jets required 

Figure 2. The number of observed events in the 5-jet bin is compared to the background expectation that is determined by 
using PYTHIA to extrapolate the number of events in data from the low jet-multiplicity control regions. The contents of the 
bins represent the number of events with 5-jets passing a given jet pt requirement. These bins are inclusive in jet pr . Results 
with various fe-tagging requirements are shown. 


Distributions for data in the inclusive > 6-jet and > 7- 
jet signal regions are shown in Fig. 4-6 compared with 
background predictions determined using extrapolations 
from three different jet-multiplicity bins. In each case, 
the distributions representing the extrapolations across 
two jet-multiplicity bins (i.e. 4 —» 6 and 5 —>■ 7) are used 
as the final background prediction whereas the other ex¬ 
trapolations are simply considered as additional valida¬ 
tion. Contributions from higher jet-multiplicity regions 


are summed to construct an inclusive sample. The sys¬ 
tematic uncertainty is constructed from the maximum 
deviation given by the various validations and for most 
signal regions is dominated by the baseline uncertainty 
obtained from the n < 5 jet regions. Results using the 
three 5-tagging selections (> 0, > 1, > 2 5-tagged jets) 
are shown in Fig. 4-6. The background systematic uncer¬ 
tainties determined from the control regions in the data 
are shown as the green shaded region in the ratio plots 
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Jet p T Requirement [GeV] 


(a) Extrapolation in 4 —> 6-jet bin 



(b) Extrapolation in 5 —> 7-jet bin with (|r/|) > 1.0 

Figure 3. The data are compared with the expected back¬ 
ground shapes in the exclusive 6- and 7-jet bins before b- 
tagging. The contents of the bins represent the number of 
events with the given number of jets passing a given jet px 
requirement. The bins with less than 10% expected signal 
contamination are control regions that are considered when 
assigning systematic uncertainties to the background yield. 
These control regions are the bins to the left of the vertical 
red lines in the plots, (a) shows the 6-jet region, and (b) 
shows the 7-jet region with (|»y|) > 1.0. 


of these figures. This procedure results in a background 
systematic uncertainty in the p!^ > 120 GeV, > 7-jet 
region of 14%, 15%, and 40% for > 0, > 1, > 2 6-tagged 
jets, respectively. 
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(a) > 6 jets, > 0 Fl agged jets 



Figure 4. The number of observed events in the inclusive 
> 6-jet (top) and > 7-jet (bottom) signal regions compared 
with expectations using the PYTHIA extrapolations from low 
jet-multiplicity control regions, as a function of the jet pt re¬ 
quirement. The distributions representing the extrapolations 
across two units in jet multiplicity (red triangles) are used as 
the final background prediction in each case, while the other 
extrapolations are treated as cross-checks. > 0 6-tagged jets 
are required. In the ratio plots the green shaded regions rep¬ 
resent the background systematic uncertainties. 


The bins in these distributions that were not as¬ 
signed as control regions represent possible signal regions, 
which may be chosen as a signal region for a particu¬ 
lar model under the optimization procedure described in 
Sec. VIIIB. The level of disagreement between the expec¬ 
tation and data is shown in Fig. 7 for the > 0 6-tagged 
jets control and signal regions. In the 6-tagged signal re¬ 
gions similar agreement is observed between data and the 
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60 80 100 120 140 160 180 200 220 240 

Jet p t Requirement [GeV] 


(a) > 6 jets, > 1 6-tagged jet 



Figure 5. Distributions shown here are as in Fig. 4 but with 
> 1 6-tagged jets required. 


predicted background, within the assigned uncertainties. 
In practice, it is seen that for most signals, the > 7-jet 
bin is preferred by the optimization procedure as a signal 
region. The data in each distribution show good agree¬ 
ment with background predictions within uncertainties. 

Systematic uncertainties on the jet-counting back¬ 
ground estimation using the extrapolation method are 
determined directly from the data as part of the back¬ 
ground validation and, by design, account for all uncer¬ 
tainties on the technique and on the reference model used 
in the projection. In contrast, systematic uncertainties 
on the signal predictions are determined from several 
sources of modeling uncertainties. The largest system¬ 
atic uncertainties are those on the background yield, the 
jet energy scale uncertainties on the signal yield (10-20% 
for most signal regions), and the uncertainty in 6-tagging 




Figure 6. Distributions shown here are as in Fig. 4 but with 
> 2 6-tagged jets required. 


efficiencies for many signal regions that require the pres¬ 
ence of 6-tagged jets (between 15-20% for signal regions 
requiring at least two 6-tags). 

An additional systematic uncertainty is included in 
these estimates in order to cover possible contamination 
of signal in the control regions for the extrapolation. The 
analysis is repeated with signal injected into the control 
regions and the backgrounds are re-computed. The re¬ 
sulting bias depends on the signal model and is found to 
be less than 5% in all cases. 

Given the good agreement between the data and the 
predictions from the jet-counting background estimation, 
there is no evidence of new physics. 
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Figure 7. Comparisons of the deviation between data and expectations in the control and signal regions without 6-tagging 
requirements are shown, as a function of the jet px requirement. The solid black line shows the relative difference between the 
observed data and the predicted background. The coarsely dashed blue distribution shows the relative systematic uncertainty 
on the background estimation. The finely dashed red distribution shows the total uncertainty on the comparison between 
background and data, including the background systematic uncertainty and all sources of statistical uncertainty from the data 
and simulation. 


VII. TOTAL-JET-MASS ANALYSIS 

A. Method and techniques 

The total-jet-mass analysis uses a topological observ¬ 
able Mf as the primary distinguishing characteristic be¬ 
tween signal and background. The observable Mf [67- 
69] is defined as the scalar sum of the masses of the four 
leading large -R jets reconstructed with a radius parame¬ 
ter R = 1.0, pt > 100 GeV and |?y| < 2.5, 

Mf= ™ jet - ( 4 ) 

P T >100 GeV 
M<2.5 

This observable was used for the first time in the -y/s = 
8 TeV search by the ATLAS Collaboration for events 


with many jets and missing transverse momentum [70] 
and provides significant sensitivity for very high-mass 
gluinos. Four-jet (or more) events are used as four large- 
R jets cover a significant portion of the central region of 
the calorimeter, and are very likely to capture most signal 
quarks within their area. This analysis focuses primarily 
on the ten-quark models mentioned in Sec. I. 

Simulation studies show that Mf provides greater sen¬ 
sitivity than variables such as Ht, the scalar sum of jet 
Pt- the masses contain angular information about the 
events by definition, whereas a variable like Ht simply 
describes the energy (or transverse momentum) in the 
event. A large Mf implies not only high energy, but also 
rich angular structure. Previous studies at the Monte 
Carlo event generator level have demonstrated the power 
of the Mf variable in the high-multiplicity events that 
this analysis targets [67, 68]. 

Figure 8 presents examples of the discrimination that 
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(a) MV for fixed m ~o = 175 GeV 

v J Xi 



Figure 8 . Comparison between signal and background for (a) 
the scalar sum of the masses of the four leading large -R jets 
Mf and (b) the difference in pseudorapidity between the two 
leading large -R jets | A 77 I. Several typical signal points are 
shown, as well as the distributions obtained from the data. 
All distributions are normalized to the same area. The selec¬ 
tion requires four or more jets, similar to the 4j regions but 
inclusive in | A 77 I. 


the Mf observable provides between the background 
(represented here by SHERPA multi-jet MC simulation) 
and several signal samples, as well as the comparison of 
the data to the SHERPA multi-jet background. Three sig¬ 


nal samples, each with to- 0 = 175 GeV and several gluino 

Xi 

masses rrig in the range 0.6 — 1.4 TeV are shown. In each 
case, the discrimination in the very high Mf region is 
similar and is dictated primarily by the gluino mass, but 
is also sensitive to the mass splitting, nig — 0 . Larger 

nig results in larger Mf, as expected. However, for the 
same mg, Mf is largest for m^o « mg/2. This is due 
to the partitioning of the energy in the final state. For 
very large to- 0 , with m-o < mg, the two quarks from 
the decay of the g are very soft and the partons from 
the decay of the Xi are relatively isotropic, slightly re¬ 
ducing the efficacy of the approach. For very low to- 0 , 
to- 0 < mg, the opposite occurs: the two quarks from 
the gluino decay have very high px and the neutralino is 
Lorentz-boosted, often to the point that the decay prod¬ 
ucts merge completely, no longer overlapping with quarks 
from other parts of the event, and the mass of the jet 
is substantially reduced . 3 In both cases, although the 
sensitivity of Mf is reduced, the overall approach still 
maintains good sensitivity. 

Another discriminating variable that is independent of 
Mf is necessary in order to define suitable control re¬ 
gions for the analysis. As in the jet-counting analysis, 
the signal is characterized by a considerably higher rate 
of central jet events as compared to the primary multi-jet 
background. This is expected due to the difference in the 
production processes that is predominantly s-channel for 
the signal, while the background can also be produced 
through u- and t-channel processes. Figure 8 addition¬ 
ally shows the distribution of the pseudorapidity differ¬ 
ence between the two leading large-1? jets, |Ap|. The 
discrimination between the signal samples and the back¬ 
ground is not nearly as significant for Ar/ as for Mf. 
However, the lack of significant correlation (Pearson lin¬ 
ear correlation coefficient of approximately 1 %) between 
the two observables makes |A?/| effective as a means to 
define additional control regions in the analysis. It is also 
observed that the shape of the distribution is relatively 
independent of the g and Xi masses and mass splittings. 

The ability of several other observables to discriminate 
between signal and background was also tested. In par¬ 
ticular, the possibility of using more detailed information 
about the substructure of jets (e.g. the subjet multiplic¬ 
ity or observables such as V-subjettiness, T 32 [71, 72]) was 
investigated. Although some additional discrimination is 
possible using more observables, these significantly com¬ 
plicated the background estimation techniques and only 
marginally increase the sensitivity of the analysis. 


O ~ Q 

^While the complete merging of the decay products of a Xi 
into a single jet may suggest that the most effective variable at low 
m-o might be the jet mass itself, typically only the lightest X? have 
enough px to be strongly collimated. Such jets thereby have very 
low jet masses. These low jet masses similar to what is expected 
from QCD radiation, making discrimination very difficult, and so 
the nominal total-jet-mass technique is maintained even in these 
regions. 
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The use of Mf in this analysis provides significant sen¬ 
sitivity as well as the opportunity to complement the 
jet-counting analysis described in Sec. VI with a fully 
data-driven background estimation that does not require 
any input from MC simulation. A template method is 
adopted in which an expected Mf distribution is con¬ 
structed using individual jet mass templates. Single-jet 
mass templates are extracted jet-by-jet from a signal- 
depleted 3-jet control region (3jCR), or training sample. 
These jet mass templates are binned in jet px and 77 , 
which effectively provides a probability density function 
that describes the relative probability for a jet with a 
given px and 77 to have a certain mass. This template 
is randomly sampled 2500 times for a single jet px and 
77 , and a precise predicted distribution of possible masses 
for the given jet is formed . 4 For an event with multiple 
jets, the jet mass templates are applied to each jet and 
the resulting predicted mass distributions are combined 
to predict the total-jet-mass Mf for that ensemble of 
jets. 

Jet mass templates are applied to jets in events in or¬ 
thogonal regions, typically with at least four large -R jets 
- the control (4jCR), validation (4jVR), and signal re¬ 
gions (4jSR) - but also in the 3jCR to test the method. 
Samples used in this way are referred to as the kinematic 
samples. The only information used is the jet px and 
77 , which are provided as inputs to the templates. The 
result is referred to as a dressed sample , which provides 
an SM prediction of the individual jet mass distributions 
for the jets in the kinematic sample. An SM prediction 
for the total-jet-mass can then be formed by combining 
the individual dressed jet mass distributions. The nor¬ 
malization of the Mf prediction - the dressed sample - 
is preserved such that the total expected yield is equal 
to the number of events in the kinematic sample. The 
procedure can be summarized as [69]: 

1. Define a control region to obtain the training sam¬ 
ple from which jet mass templates are to be con¬ 
structed; 

2. Derive a jet mass template binned in jet 77 and px 
using a smoothed Gaussian kernel technique; 

3. Define a kinematic sample as either another control 
region or the signal region; 

4. Convolve the jet mass template with the kinematic 
sample using only the jet px and 77 ; 

5. Obtain a sample of dressed events which provides 
the data-driven background estimate of Mf. 

The key assumption in this approach is that the jet 
kinematics factorize and are independent of the other jets 


4 2500 times was found to be the best balance between the pre¬ 
cision of the result and computational time. 


in the event. Deviations from this approximation may oc¬ 
cur due to effects that are not included in the derivation 
of the jet mass templates. In particular, the composi¬ 
tion of quarks and gluons can vary across different sam¬ 
ples [73], and quark and gluon jets have been observed 
to have different radial energy distributions [74]. Other 
experimental affects, arising from close-by or overlapping 
jets, can also have an effect. For this reason, extensive 
tests are performed in the 4jCR and 4jVR, as defined 
in Sec. VIIB, to estimate the size of the correction fac¬ 
tors needed to account for any sample dependence, and 
to assess systematic uncertainties. The entire procedure 
is tested first in SHERPA multi-jet MC simulation, which 
shows minimal differences between the template predic¬ 
tion and observed mass spectrum. 


Region 

77 jet 

|At?| 

Px 

Pt 

Mf 

Name 



[GeV] 

[GeV] 

[GeV] 

3jCR 

CO 

II 

<U 

- 


- 

- 

4jCR 

Al 

CD 

> 1.40 

> 100 

> 250 

> 100 

: 

4jVR 

^-jet ^ 4 

1.0-1.40 

> 100 

> 250 

> 100 

: 

SRI 



> 250 


> 625 

SR100 

TT-jet ^ 4 

< 0.7 

> 100 

> 100 

> 350 (binned) 

SR.250 



> 250 


> 350 (binned) 


Table I. Control (CR), validation (VR), and signal regions 
(SR) used for the analysis. pf and p% represent the transverse 
momentum of the third and fourth jet in px, respectively. 


B. Signal and control region definitions 

The Mf and |A? 7 | observables form the basis for the 
signal region definition for the analysis, where IA 77 I is 
used to define control regions for testing the background 
estimation in data. A requirement of | A 77 I < 0.7 is found 
to have the best signal sensitivity over the entire plane 
of (my, m^o). In this optimization, the background con¬ 
tribution is modeled by multi-jet events simulated with 
SHERPA. 

An optimization study indicated that when using a 
single Mf selection, Mf > 625 GeV provides the best 
sensitivity to many signal hypotheses, and gives the best 
expected sensitivity at high my. A single-bin signal re¬ 
gion (SRI) is therefore defined with Mf > 625 GeV and 
a 250 GeV px threshold applied to the third leading in 
px large-i? jet. This region has an acceptance of 0.26% 
for the my = 600 GeV, tti-o = 50 GeV signal point. This 
acceptance grows rapidly with gluino mass to 11 % for 
the point my = 1000 GeV, 777.-0 = 600 GeV, and is only 
weakly dependent on the neutralino mass. 
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Figure 9. (a) Total-jet-mass in the 4jVR with px > 100 GeV. 
The reweighted template is shown in the hatched blue his¬ 
togram. (b) Total-jet-mass in the 4jVR with Px > 250 GeV. 
The 4jVR Mf spectra are shown in the open black squares. 
The total systematic uncertainty due to the smoothing pro¬ 
cedure, finite statistics in control regions, and the difference 
between template prediction and the data observed in the 
4jCR is shown in green. 


Figure 10. Total-jet-mass in the 4jSR (a) using pf > 100 GeV 
(SR100) and (b) using px > 250 GeV (SR250). For the SR100 
selection, the reweighted template (built in the 3jCR, and 
reweighted jet-by-jet in the 4jCR) is shown in the hatched 
blue histogram. The total systematic uncertainty due to the 
smoothing procedure, finite statistics in control regions, and 
the difference between template prediction and the data ob¬ 
served in the 4jCR is shown in green. 


A second set of signal regions is used to further improve 
the power of the analysis by making use of the shape of 
the Mf distribution. Two selections on the third lead¬ 
ing jet in px (pt) are use d, pf > 100 GeV (SR100) and 
pf > 250 GeV (SR250). This provides better sensitiv¬ 
ity to the full range of gluino masses considered, com¬ 
pared to SRI. The lower px region, SR100, has better 
sensitivity for lower gluino masses, whereas SR250 has 
improved sensitivity for higher masses. All other selec¬ 
tions are unchanged. In this case, a lower threshold of 


Mf > 350 GeV is used and the observed data are com¬ 
pared to the template predictions in bins of Mf. The 
improvements in the sensitivity obtained by adding these 
additional signal regions and using the shape of the Mf 
spectrum are described below. The full set of selection 
criteria is listed in Tab. I. 

The jet multiplicity and |Ap| are used to define the 
control regions. The 3jCR, with exactly three jets, is 
used to train the background templates previously dis¬ 
cussed. In the remaining control and validation regions, 
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each requiring > 4 jets, the |Ary| selection suppresses the 
signal contribution and is used to define the 4jCR and 
4jVR. In the > 4-jet regions, the |A 77 I selection value for 
the control regions is chosen to be larger than an inver¬ 
sion of the signal region selection, resulting in the selec¬ 
tions presented in Tab. I. These control region definitions 
permit studies of the full Mf spectrum as well as com¬ 
parisons of data and SM predictions without significant 
signal contamination. 

C. Validation and systematic uncertainties 

Many tests are performed using the 3jCR as both the 
training sample and the kinematic sample in order to 
determine the robustness of the method. The selection 
requires that there be exactly three large- 1 ? jets in the 
event, as described in Tab. I. The dependence of the tem¬ 
plate on the jet in question (leading, subleading, etc) is 
tested, as well as the dependence of the template on the 
jet kinematics. It is determined that it is optimal to 
define separate templates for each of the three jet cate¬ 
gories (leading, subleading, and third jet) and to bin the 
templates according to the jet pr and 77. 5 In the 4-jet 
regions, the fourth jet uses the template derived for the 
third jet in the 3jCR: tests in the 4jCR and 4jVR in¬ 
dicate very good agreement between this template and 
the observed spectrum. As a first test, the Mf template 
constructed from the 3-jet kinematic sample is compared 
to the actual Mf distribution in 3-jet events, and very 
good agreement is observed. 

There are two intrinsic sources of systematic uncer¬ 
tainty associated with the template procedure: the uncer¬ 
tainty due to finite statistics in the 3jCR training sample 
(the variance), and the uncertainty due to the smooth¬ 
ing procedure in the template derivation (the bias). The 
former is estimated by generating an ensemble of Mf 
templates and taking the ± 1 <t deviations (defined as the 
±34% quantile) with respect to the median of those vari¬ 
ations as the uncertainty, bin by bin. The systematic un¬ 
certainty due to the smoothing procedure is determined 
using the fact that a Gaussian kernel smoothing is ap¬ 
plied to the template. The full difference between the 
nominal template and a template constructed using a 
leading-order correction for the bias, derived analytically 
in Ref. [69], is taken as the systematic uncertainty. The 
systematic uncertainty due to finite control region statis¬ 
tics is chosen to be larger (by setting the size of the kernel 
smoothing) than that due to the smoothing procedure 
since the former is more accurately estimated. 

A small level of disagreement (between 5 to 15%) is 
observed when comparing the observed mass to the pre¬ 
dicted mass in the 4jCR: a reweighting derived in the 


5 It is observed that the difference between the leading and sub¬ 
leading jet templates is minimal, but that the third jet exhibits 
qualitatively different masses as function of the jet px- 


4jCR (as a function of each individual jet mass) is then 
applied to the individual jet masses prior to the construc¬ 
tion of the Mf for each event. After the reweighting 
the agreement is substantially improved at high total- 
jet-mass. Figure 9 presents the total-jet-mass Mf in 
the 4jVR using pf > 100 GeV. The reweighted template 
agrees very well with the observed Mf distribution in the 
4jVR— a sample completely independent from where the 
reweighting was derived— validating both the template 
method and the reweighting. The full magnitude of the 
reweighting on the total-jet-mass distribution is taken as 
a systematic uncertainty of the method. The total sys¬ 
tematic on the background prediction therefore includes 
both the intrinsic systematic uncertainty given by the 
variance and the bias, as well as the difference observed 
in the 4jCR. The Mf distribution is also shown for the 
4j VR for the case in which pf > 250 GeV. No reweighting 
is required when using the significantly higher pf selec¬ 
tion since the observed effects due to topological differ¬ 
ences in the training sample compared to the kinematic 
sample are suppressed. In order to account for any re¬ 
maining disagreement, the difference between the data 
and template prediction in the 4jCR is applied as a fur¬ 
ther systematic. The total uncertainty therefore includes 
again both the instrinsic background estimation uncer¬ 
tainties and the disagreement observed in the 4jCR. 

One possible concern for the template technique is that 
it assumes that the same mechanism is responsible for 
generating the individual jet masses in both the con¬ 
trol and signal regions. In order to test the extent to 
which a different composition of processes may affect the 
derived templates, the assumption that multi-jet events 
are the only background in the 3jCR and 4j regions is 
modified by injecting separately a sample of SHERPA tt 
MC simulation events (assuming SM cross-sections) into 
the full procedure. The resulting background estimates 
are fully consistent with the prediction without the in¬ 
jection - indicating that the technique is not sensitive 
to contamination from top quark production - and thus 
no additional systematic uncertainty is assessed for the 
potential presence of specific background processes. A 
similar procedure is performed for signal processes (as¬ 
suming standard g production cross-sections) and again 
no impact of signal contamination on the constructed 
background templates is observed. 

Figure 10 shows the total-jet-mass in the 4jSR com¬ 
pared to the template prediction. For both SR100 
and SR250, the total systematic error on the template 
method is also shown in the ratio plot in the lower panel 
of each distribution. The template predictions are clearly 
consistent with the observed data. Thus there is no in¬ 
dication of new physics in these results. 

Systematic uncertainties associated with the scale and 
resolution of large-1? jet mass and energy [65] are signif¬ 
icantly reduced by the use of a data-driven background 
estimate: residual effects may remain due to differences 
between the 3jCR and the 4j regions, and these are re¬ 
flected in the systematic uncertainties assessed by the 
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difference between the template prediction and observed 
spectrum in the 4jCR. The uncertainties due to the back¬ 
ground estimation method are dominated by propaga¬ 
tion of the statistical uncertainty from the 3jCR: these 
are typically 5-10%, except in the highest Mf bins of 
SR100 and SR250, where they can extend to 20-40%. In 
addition, the observed difference systematic uncertainty 
from the 4jCR varies from 5% to 15%. Signal reconstruc¬ 
tion - both in terms of selection efficiency and the Mf 
spectrum predicted for a given m.g , o combination - is 
sensitive to the kinematic uncertainties associated with 
the final state jets in the analysis. The impacts of these 
systematic uncertainties are directly assessed by varying 
the kinematics within the uncertainties and reported in 
Sec. VIII. Jet mass scale uncertainties have the largest 
effect, which for SRI range from 30% for very low rrig to 
15% for very high rrig. In the cases of SR100 and SR250, 
the impact of the jet mass scale uncertainty also domi¬ 
nates, and varies across the Mf spectrum from 10-20% 
at lower Mf up to 50% for the very highest Mf bin in 
the spectrum for low rrig. The luminosity uncertainty of 
3% also affects the signal only. 


VIII. RESULTS AND INTERPRETATIONS 

As no significant excess is observed in data in either 
analysis, a procedure to set limits on the models of in¬ 
terest is performed. A profile likelihood ratio combining 
Poisson probabilities for signal and background is com¬ 
puted to determine the confidence level (CL) for consis¬ 
tency of the data with the signal-plus-background hy¬ 
pothesis (CL s+ b). A similar calculation is performed for 
the background-only hypothesis (CLf,). From the ratio 
of these two quantities, the confidence level for the pres¬ 
ence of signal (CL S ) is determined [75]. Systematic un¬ 
certainties are treated via nuisance parameters assum¬ 
ing Gaussian distributions. In all cases, the nominal 
signal cross-section and uncertainty are taken from an 
envelope of cross-section predictions using different PDF 
sets and factorization and renormalization scales, as de¬ 
scribed in Ref. [76]. As discussed in Sec. IV, the region 
with ( nig — o) < 100 GeV is not considered in this 
analysis in order to ensure that the results are insensi¬ 
tive to the effects of ISR, since the uncertainties cannot 
be assessed for the UDD decays considered here. 

The total-jet-mass analysis is designed to be agnostic 
to the flavor composition of the signal process and to re¬ 
move any reliance on MC simulations of these complex 
hadronic final states. The jet-counting analysis provides 
the opportunity to enhance sensitivity to specific heavy- 
flavor compositions in the final state and to explore var¬ 
ious assumptions on the branching ratios of the bench¬ 
mark signal processes studied in this paper. The results 
obtained from the total-jet-mass analysis in the inclusive 
final state are presented first, and then the specific sen¬ 
sitivity provided by the jet-counting analysis to the full 
branching ratio space is presented. 
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Figure 11. Expected and observed exclusion limits in the 
(rrig, o) plane for the ten-quark model given by the total- 
jet-mass analysis. Limits are obtained by using the signal 
region with the best expected sensitivity at each point. The 
dashed black lines show the expected limits at 95% CL, with 
the light (yellow) bands indicating the la excursions due to 
experimental and background-only theory uncertainties. Ob¬ 
served limits are indicated by medium dark (maroon) curves, 
where the solid contour represents the nominal limit, and the 
dotted lines are obtained by varying the signal cross-section 
by the renormalization and factorization scale and PDF un¬ 
certainties. 


A. Total-jet-mass analysis 

The observed and expected event yields are presented 
in Tab. II, III, and IV for the three signal regions SRI, 
SR100 and SR250 respectively. The single-bin signal re¬ 
gion selection (SRI) is reported in addition to the binned 
Mf results in SR100 and SR250 in order to provide yields 
that can be easily reinterpreted for other signal hypothe¬ 
ses. In the case of the binned Mf signal regions, a binned 
fit (where the number and size of the bins were optimized) 
is performed that takes into account the predictions for 
each Mf range. This approach provides greater sensi¬ 
tivity to small deviations from the template predictions. 
The correlation of the uncertainties in the bins of the 
Mf spectrum are accounted for by evaluating the full 
correlation matrix. The result leads the analysis to treat 
the different bins as fully uncorrelated for the variance, 
which is the largest component of the background uncer¬ 
tainties. All other uncertainties treat the bins of the Mf 
spectrum as fully correlated. 

Figure 11 shows both the expected and observed 95% 
CL limits in the (rrig, m^o) mass plane when the signal 
region that provides the best expected exclusion is used 
for each mass combination. The dashed black line shows 
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Summary yield table for SRI 


Mf Bin 

Expected SM 

Obs. 

m,g = 600 GeV 
m -o = 50 GeV 

?rig = 1 TeV 
m~ o = 600 GeV 

m.g = 1.4 TeV 

m- o = 900 GeV 

> 625 GeV 

160±9.7 ±|2 

176 

70±4.2 ±25±30 (0.26%) 

55T0.51 ±8.6 ±14 (11%) 

6.3±0.07 ±0.46±2.5 (35%) 


Table II. Table showing the predicted in the SM and observed number of events in SRI as well as three representative signal 
scenarios. Acceptances (including efficiency) of the various signals are listed in parentheses. The background uncertainties are 
displayed as statistical + systematic; the signal uncertainties are displayed as statistical ± systematic ± theoretical. 


Summary yield table for SR100 


Mf Bin 

Expected SM 

Obs. 

irig = 600 GeV 
m~o = 50 GeV 

rrig = 1 TeV 

m- o = 600 GeV 

rrig = 1.4 TeV 
m~ o = 900 GeV 

350 - 400 GeV 

4300±78 

5034 

200±7.2±22±35 

5.8±0.17±1.3±1.5 

0.19±0.01±0.04±0.07 

400 - 450 GeV 

2600±49 ±||g 

2474 

200±7.1±9.5±35 

9.7±0.21±2.2±2.5 

0.31±0.02±0.07±0.12 

450 - 525 GeV 

2100±42 +Ho 

1844 

280±8.4±13±49 

26±0.35±4.3±6.7 

0.88±0.03±0.14±.34 

525 - 725 GeV 

960±25 t 2 2 ° 0 ° 0 

1070 

280±8.4±57±49 

77±0.60±3.2 

3.6±0.05±0.36±1.4 

> 725 GeV 

71±7.0 tj? 

79 

35.±2.9±18±6.0 

35±0.40±9.9±9.0 

4.8±0.06±0.61±1.9 


Table III. Table showing the predicted in the SM and observed number of events in SR100 as well as three representative signal 
scenarios. The background uncertainties are displayed as statistical + systematic; the signal uncertainties are displayed as 
statistical ± systematic ± theoretical. 


the expected exclusion limits, and the yellow band rep¬ 
resents the experimental uncertainties on this limit. The 
solid line shows the observed limit, with the finely dashed 
lines indicating the ±ler variations due to theoretical un¬ 
certainties on the signal production cross-section given 
by renormalization and factorization scale and PDF un¬ 
certainties. All mass limits are reported conservatively 
assuming the — l^fhemry sig na l production cross-section. 
At low m^o, the region with gluino mass rrig < 750 GeV 
is excluded. Excluded nig masses rise with increas¬ 
ing m- o, up to a maximum exclusion of approximately 
TYin < 870 GeV at m~o = 600 GeV. No models with 

y Xi 

m < o > 650 GeV are excluded. 

X1 

B. Jet-counting analysis 

In order to set limits on individual branching ratios, 
it is necessary to refer to the structure of the couplings 
that are allowed. From Eq. (1), it is clear that each 
RPV decay produces exactly two down-type quarks of 
different flavor and one up-type quark. Since the cross- 
section for gluino production is not dependent upon the 
X” ik parameters, it is not possible to directly probe or 
set limits upon any individual A"- fc parameter. Instead, 
results are categorized based upon the probability for an 
RPV decay to produce a t-quark, a 5-quark, or a c-quark. 
These branching ratios are denoted by BR(f), BR(6), and 
BR(c), respectively. These branching ratios are partially 
constrained. The branching ratios for decays including 
u-, c-, and t-quarks (all given by the flavor index i in the 
A" fc couplings) must sum to one and must therefore sat¬ 


isfy BR(f) + BR(c) < 1. The branching ratios to decays 
including each down-type quark (as given by the flavor 
indices j and k in the A"- fe couplings) are independent 
of the up-type branching ratios. At most, one 5-quark 
can be produced in such an RPV decay. Simultaneous 
nonzero A” values can result in nontrivial branching ra¬ 
tio combinations. 

Results using the jet-counting analysis are determined 
for different hypotheses on the branching ratios of RPV 
decays to t, 5, c, and light-flavor quarks. The selec¬ 
tion requirements for the signal regions are optimized 
separately for each of these hypotheses. When running 
the optimization, the full limit-setting procedure is per¬ 
formed under the assumption that the expected number 
of background events is observed in the data, taking all 
statistical and systematic uncertainties into account. The 
results of this optimization are provided in Tab. V. The 
first portion of Tab. V shows the optimization results 
and the comparison of the data with background predic¬ 
tions for the six-quark signal models under the assump¬ 
tion that (BR(i), BR(5), BR(c))=(0%, 0%, 0%). In this 
simple model, it is equivalent to say that only the term 
given by A , / 12 is nonzero. Explicitly, this flavor hypothe¬ 
sis forces the RPV decays to result only in light quarks. 
Below this, the table shows the same comparisons under 
the assumption that (BR(f), BR(5), BR(c))=(0%, 100%, 
0%) corresponding to only RPV terms given by A" 13 and 
Ai 23 - The second half of Tab. V is analogous to the first, 
only with BR(t)=100%. The signal acceptance is largely 
affected by BR(t) and BR(5) due to the presence of signal 
regions with 5-tagged jets. Because this search requires 
many high-px jets, increased BR(t) results in a lower 























17 


Summary yield table for SR250 


Mf Bin 

Expected SM 

Obs. 

rrig = 600 GeV 
m~ o =50 GeV 

rrig = 1 TeV 
m~ o = 600 GeV 

nig = 1.4 TeV 

m- o = 900 GeV 

350 - 400 GeV 

1400T35 

1543 

83±4.6 ±15±14 

3.3±0.12 ±0.78±0.85 

0.17±0.01 ±0.03±0.07 

400 - 450 GeV 

920±33 

980 

92±4.8 ±11±16 

5.6±0.16 ±1.5±1.5 

0.27±0.01 ±0.07±0.11 

450 - 525 GeV 

780±33 ±94 

823 

140±5.8 ±15±23 

17±0.28 ±3.3±4.4 

0.79±0.02 ±0.13±0.31 

525 - 725 GeV 

490±24 

495 

160T6.2 ±30.±27 

56±0.51 ±4.1±15 

3.3±0.05 ±0.34±1.3 

> 725 GeV 

37±5.5 

42 

22±2.3 ±9.1±3.9 

27±0.36 ±7.4±7.0 

4.4±0.06 ±0.56±1.7 


Table IV. Table showing the predicted in the SM and observed number of events in SR250 as well as three representative signal 
scenarios. The background uncertainties are displayed as statistical ± systematic; the signal uncertainties are displayed as 
statistical ± systematic + theoretical. 






Figure 12. Expected and observed cross-section limits for 
the six-quark gluino models for (a) the case where no gluinos 
decay into heavy-flavor quarks, and (b) the case where every 
gluino decays into a 6-quark in the final state. 


Figure 13. Expected and observed cross-section limits for the 
six-quark gluino models for (a) the case where each gluino is 
required to decay into a top quarks, and (b) the case where 
every gluino decays into a b-quark and a top quark. 
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pp->gg->6q, BR(c)=0% 

All limits at 95% CL i/i = 8 TeV, 20.3 fb' 1 
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Figure 14. Expected and observed mass exclusions at the 95% 
CL in the BR(f) vs. BR(fe) space for BR(c)=0%. Each point 
in this space is individually optimized and fit. Masses below 
these values are excluded in the six-quark model. Bin centers 
correspond to evaluated models. 
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All limits at 95% CL \/s = 8 TeV, 20.3 fb' 1 
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Figure 15. Expected and observed mass exclusions at the 95% 
CL in the BR(t) vs. BR(fe) space for BR(c)=50%. Each point 
in this space is individually optimized and fit. Masses below 
these values are excluded in the six-quark model. Bin centers 
correspond to evaluated models. 


acceptance from larger energy sharing in a higher multi¬ 
plicity of final state objects. For this reason, the corners 
of the BR(f) vs. BR(6) space are shown here. Since 
the sensitivity to increased BR(c) comes from ^-tagging 
configurations that are designed to efficiently select b- 
jets, the effect on the signal acceptance is dominated by 
BR(6). For this reason, the focus of this discussion is on 
the BR(6) degree of freedom. However, several results 
with various values of BR(c) are presented below. 

The results of performing the limit-setting procedure 
on the data in the signal regions are shown in Fig. 12 and 
13 for various flavor branching ratio hypotheses as a func¬ 
tion of gluino mass for the six-quark model. These results 


show both the expected and observed cross-section lim¬ 
its in comparison to the predicted cross-section from the 
theory. Under the assumption that all RPV decays are to 
light-flavor quarks (BR(6)=BR(t)=BR(c)=0%), gluino 
masses of rrig < 853 GeV (expected) and < 917 GeV 
(observed) are excluded at the 95% CL. Alternatively for 
the scenario where BR(6)=100% while the other heavy- 
flavor branching ratios are zero, exclusions of m.g < 921 
GeV (expected) and trig < 929 GeV (observed) are found. 
Similarly, for the case where BR(6)=BR(f)=100%, exclu¬ 
sions of uig < 938 GeV (expected) and trig < 874 GeV 
(observed) are found. More generally, excluded masses 
as a function of the branching ratios of the decays are 
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g-g production, g-> qq^, qqq 
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m g [GeV] 



Sample 

mg 

[GeV] 

Jet pt 

req. 

[GeV] 

# of 
jets 

# of 
6 -tags 

Signal 

(Acceptance) 

Back¬ 

ground 

Data 



(BR(t), BR(6), BR(c))=(0%, 0%, 0%) 


500 

120 

7 

0 

600=4230 (0.7%) 

370=460 

444 


600 

120 

7 

0 

410=4100 (1.5%) 

370=460 

444 


800 

180 

7 

0 

13=44 (0.4%) 

6.1=42.2 

4 


1000 

180 

7 

0 

6.8=42.3 (1.4%) 

6.1=42.2 

4 


1200 

180 

7 

0 

2.7=40.5 (3.0%) 

6.1=42.2 

4 


(BR(t), BR(6), BR(c))=(0%, 100%, 0%) 


500 

80 

7 

2 

1900=4400 (2.1%) 

1670=4190 

1560 


600 

120 

7 

1 

300=460 (1.1%) 

138=426 

178 


800 

120 

7 

1 

131±25 (4.1%) 

138=426 

178 


1000 

180 

7 

1 

4.4=41.0 (0.9%) 

2.3=41.0 

1 


1200 

180 

7 

1 

1.86=40.31 (2.1%) 

2.3=41.0 

1 

I4UU 

iVl 

(BR(i), BR(b), BR(c))=(100%, 0%, 0%) 



500 

80 

7 

1 

4600T800 (5.0%) 

5900=4700 

5800 


600 

100 

7 

1 

940=4190 (3.5%) 

940=4140 

936 


800 

120 

7 

1 

1084=18 (3.4%) 

138=426 

178 


1000 

120 

7 

1 

42=46 (8.5%) 

138=426 

178 


1200 

180 

7 

1 

1.3±0.4 (1.5%) 

2.3=41.0 

1 


(BR(t), BR(6), BR(c)) = (100%, 100%, 0%) 


500 

80 

7 

2 

3600=4600 (3.9%) 

1670=4190 

1560 


600 

80 

7 

2 

2300=4400 (8.6%) 

1670=4190 

1560 


800 

120 

7 

2 

94=415 (3.0%) 

38=417 

56 


1000 

120 

7 

2 

37=46 (7.5%) 

38=417 

56 


1200 

140 

7 

2 

5.5=41.0 (6.2%) 

10=45 
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Table V. Requirements as optimized for the six-quark model 
under a variety of gluino mass hypotheses when the RPV ver¬ 
tex has various branching ratio combinations corresponding 
to respective RPV terms given by A” k being nonzero. The op¬ 
timized signal region selection requirements are shown along 
with the resulting background and signal expectations and 
the number of observed data events. The nominal signal ac¬ 
ceptance (including efficiency) is also shown for each result. 
Quoted errors represent both the statistical and systematic 
uncertainties added in quadrature. 


(b) With 6-tagging optimization 

Figure 16. Expected and observed exclusion limits in the 
( nig , m^o) plane for the ten-quark model given by the jet¬ 
counting analysis, (a) shows the results when the branching 
ratios for the RPV decay are considered inclusively, without 
any 6-tagging requirements applied. This figure is analogous 
to Fig. 11. (b) shows the exclusion results when b-tagging 
requirements are allowed to enter into the optimization pro¬ 
cedure, improving limits significantly. 


presented in Fig. 14 and 15 where each bin shows the 
maximum gluino mass that is excluded for the given de¬ 
cay mode. 

The event selection is optimized separately for the ten- 
quark model. Table VI shows the results for the ten- 
quark model with all UDD couplings allowed, as in the 


total-jet-mass analysis, when the number of 6-tagged jets 
is also used as a variable in the optimization procedure. 
For the flavor-agnostic model where all couplings are 
equal, Fig. 16(a) shows both the expected and observed 
limits in the (mg, m^o) mass plane when the signal re¬ 
gion that provides the best expected exclusion is used 
for each mass point, not including signal regions contain¬ 
ing 6-tagged jets. The shapes of the contours are given 
by discontinuous changes in the optimized signal regions 
and fluctuations well within the given uncertainties. At 
m-o ~ 100 GeV, models with mg < 700 GeV are ex¬ 
cluded. Figure 16(b) shows the exclusion when signal re¬ 
gions with 6-tagged jets are considered as part of the op¬ 
timization and increase the sensitivity up to mg < 1 TeV 
for moderate mg — m^o mass splittings. In the ten-quark 
model, there is a significant probability that the cascade 
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Sample 
(mj, m-o) 

Jet pt req. 
[GeV] 

# jets 

# 6-tagged jets 

Signal 

(Acceptance) 

Background 

Data 

(400 GeV, 50 GeV) 

80 

7 

2 

1900T400 (0.5%) 

1670T190 

1558 

(400 GeV, 300 GeV) 

80 

7 

2 

2500T600 (0.7%) 

1670T290 

1558 

(600 GeV, 50 GeV) 

120 

7 

1 

180±40 (0.7%) 

138T26 

178 

(600 GeV, 300 GeV) 

80 

7 

2 

2200T350 (8.3%) 

1670T200 

1558 

(800 GeV, 50 GeV) 

120 

7 

1 

95±16 (3.0%) 

138T26 

178 

(800 GeV, 300 GeV) 

120 

7 

1 

172±28 (5.4%) 

138T26 

178 

(800 GeV, 600 GeV) 

120 

7 

1 

150±23 (4.7%) 

138T26 

178 

(1000 GeV, 50 GeV) 

220 

6 

1 

7.Oil.3 (1.4%) 

3.8T3.0 

5 

(1000 GeV, 300 GeV) 

120 

7 

1 

67±8 (14%) 

138T26 

178 

(1000 GeV, 600 GeV) 

120 

7 

1 

101T13 (20%) 

138T26 

178 

(1000 GeV, 900 GeV) 

120 

7 

1 

33±4 (6.7%) 

138A26 

178 

(1200 GeV, 50 GeV) 

220 

6 

1 

3.8T0.7 (4.3%) 

3.8T3.0 

5 

(1200 GeV, 300 GeV) 

180 

7 

1 

2.01T0.32 (2.3%) 

2.3T1.0 

1 

(1200 GeV, 600 GeV) 

140 

7 

1 

18.9T2.3 (21%) 

41±12 

45 

(1200 GeV, 900 GeV) 

140 

7 

1 

12.6T1.5 (14%) 

41±12 

45 


Table VI. Requirements as optimized for the ten-quark model under a variety of mass hypotheses when all couplings are 
nonzero and equal and 6-tagging requirements are considered as part of the optimization procedure. The optimized signal region 
selection requirements are shown along with the resulting background and signal expectations and the number of observed data 
events. The nominal signal acceptance (including efficiency) is also shown for each result. Quoted errors represent both the 
statistical and systematic uncertainties added in quadrature. 


decays of the gluinos produce at least one 6- or f-quark, 
and so the requirement of a 6-tagged jet improves the 
sensitivity of the analysis. 


C. Comparisons 

Model-independent upper limits on non-SM contribu¬ 
tions are derived separately for each analysis, using the 
SRI signal region for the total-jet-mass analysis. A set 
of generic signal models, each of which contributes only 
to the individual signal region, is assumed and no ex¬ 
perimental or theoretical signal systematic uncertainties 
are assigned other than the luminosity uncertainty. A 
fit is performed in the signal regions to determine the 
maximum number of signal events which would still be 
consistent with the background estimate. The resulting 
limits on the number of non-SM events and on the visible 
signal cross-section are shown in the rightmost columns 
of Tab. VII. The visible signal cross-section (tr v i s ) is de¬ 
fined as the product of acceptance (A), reconstruction 
efficiency (e) and production cross-section (fT pro d); it is 
obtained by dividing the upper limit on the number of 
non-SM events by the integrated luminosity. The results 
of these fits are provided in Tab. VII. 

The interpretations of the results of the jet-counting 
and total-jet-mass analyses are displayed together in 
Fig. 17 for the ten-quark model. This figure allows for the 
direct comparison of the results of the various analyses. 
Without 6-tagging requirements, the jet-counting analy- 



Figure 17. Expected and observed exclusion limits in the 
(m g, m-o) plane for the ten-quark model for the jet-counting 
analysis (with and without 6-tagged jets) and the total-jet- 
mass analysis. 


sis sets slightly lower expected limits than the total-jet- 
mass analysis. With 6-tagging requirements, the limits 
are stronger for the jet-counting analysis. The observed 
limits from the total-jet-mass analysis and jet-counting 
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Signal Region 

Expected 

Obs. 

Fo 

•^non-SM 

Exp. 

•^non-SM 

Obv. 

CT vis [fb] 
Exp. 

cVis [fb] 
Obv. 


SRI ( Mf ) 

160^3® 

176 

0.39 

49 

64 

2.4 

3.2 

(^-jetj Ft ■> ^b— tags) = 

(7, 120 GeV, 0) 

370 ± 60 

444 

0.12 

44 

38 

2.2 

1.9 

(^-jetj Ft 5 f^b— tags) = 

(7, 180 GeV, 0) 

6.1 ± 2.2 

4 

> 0.5 

19 

10 

0.9 

0.5 

(^-jetj Ft i ^b— tags) = 

(7, 120 GeV, 1) 

138 ± 26 

178 

0.09 

56 

42 

2.8 

2.1 

(^-jetj Ft ’ ^b— tags) = 

(7, 180 GeV, 1) 

2.3 ± 1.0 

1 

> 0.5 

4 

4 

0.2 

0.2 

(^-jetj Ft ’ n b— tags) = 

(7, 80 GeV, 2) 

1670 ± 190 

1560 

> 0.5 

38 

38 

1.9 

1.9 

(^-jetj Ft ’ ^b— tags) = 

(7, 120 GeV, 2) 

38 ± 17 

56 

0.17 

36 

52 

1.8 

2.6 


Table VII. Table showing upper limits on the number of events and visible cross sections in various signal regions. Columns two 
and three show the expected and observed numbers of events. The uncertainties on the expected yields represent systematic 
and statistical uncertainties. Column four shows the probabilities, represented by the po values, that the observed numbers of 
events are compatible with the background-only hypothesis (the po values are obtained with pseudo-experiments). Columns 
five and six show respectively the expected and observed 95% CL upper limit on non-SM events (Al non -SM), and columns seven 
and eight show respectively the 95% CL upper limit on the visible signal cross-section (cr v i s = a pl 0 a x A x t = iV non _sM/'£)• In 
the case where ALxpected exceeds ^observed, Po is set to > 0.5. 


analysis with 6-tagging requirements are also compara¬ 
ble. ‘ 


IX. CONCLUSIONS 

A search is presented for heavy particles decaying into 
complex multi-jet final states using an integrated lumi¬ 
nosity of 20.3 ± 0.6 fb -1 of y/s = 8 TeV pp collisions 
with the ATLAS detector at the LHC. Two strategies 
are used for both background estimation and signal dis¬ 
crimination. An inclusive data-driven analysis using the 
total-jet-mass with a template method for background 
estimation is performed as well as a jet-counting anal¬ 
ysis that includes exclusive heavy-flavor signal regions 
and provides limits on different branching ratios for the 
benchmark SUSY RPV UDD decays. For the ten-quark 
model, results from both analyses are presented with 
comparable conclusions. When the jet-counting analy¬ 
sis includes sensitivity to heavy flavor given by 6-tagging 
requirements, mass exclusions are further increased. 

Exclusion limits at the 95% CL are set extending up to 
m.g = 917 GeV in the case of pair-produced gluino decays 
to six light quarks and up to rrig = 1 TeV in the case 
of cascade decays to ten quarks for moderate rrig — m-o 
mass splittings. Limits are also set on different branching 
ratios by accounting for all possible decay modes allowed 
by the A"- fc couplings in full generality in the context of R- 
parity-violating supersymmetry. These results represent 
the first direct limits on many of the models considered 
as well as the most stringent direct limits to date on those 
models previously considered by other analyses. 
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